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Abstract 

We study the four channels associated with neutrino-deuteron breakup reac- 
tions at next-to-next to leading order in effective field theory. We find that the 
total cross-section is indeed converging for neutrino energies up to 20 MeV, 
and thus our calculations can provide constraints on theoretical uncertainties 
for the Sudbury Neutrino Observatory. We stress the importance of a direct 
experimental measurement to high precision in at least one channel, in order 
to fix an axial two-body counterterm. 
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I. INTRODUCTION 



Our understanding of electroweak processes on the deuteron is reaching a critical junc- 
ture. The Sudbury Neutrino Observatory is taking data, and a thorough understanding of 
neutrino-deuteron scattering is an important part of the analysis in that experiment |jl|,0] . Of 
further note, there is a proposal for a high-precision measurement of the reaction Ued e~pp 
in the ORLAND experiment @]. 

Theorists have made a tremendous effort to understand i'{v)-d scattering in a potential 
model framework fl-[lO|, with the most recent independent efforts agreeing within a few 
percent at low-energy JTT|JT^ . Given the ongoing experimental interest in these processes, 
efforts began to study vij7)-d scattering in the language of effective field theory (EFT) 
for neutrino energies below 20 MeV. The EFT work employed the power-counting scheme 
of Kaplan, Savage, and Wise ||14| and included pions. In working to next-to-leading order 
(NLO), it was found in ref. 



vd — > 

LlA- 



that theoretical uncertainties in the reactions vd — > i/np, 
Vnp, and V^d e^nn were dominated by an unknown axial two-body counterterm 
It was possible to find different values of Li^a which provided excellent fits to the 

Further, we found that the ratio of charged- 



potential model calculations of refs. |]10|,|TT 



current to neutral current (CC/NC) cross-sections was insensitive to this counterterm. This 
confirmed the insensitivity to short- distance physics first discussed by ref. [11]. However, 



it is not clear whether a power-counting scheme exists which would allow us to extend the 
theory with pions to higher-order p!5Hl8|, as would be required to constrain our theoretical 
uncertainties. 

In this work, we employ the theory without pions |jl9| (also see earlier works |ll^ , ^0| -p6| ) 
which has proven so successful in providing for high-precision calculations for other processes 
like np — (i7 [^^^. This allows us to consider next-to-next-to leading order (NNLO) cor- 



rections, and to incorporate Coulomb effects cleanly, following the prescription of Kong and 
Ravndal ^9] with some generalization. As a result, we present a complete set of calculations 
for all four reaction channels including the most-important pp channel. No new parameters 
are introduced at this order, so it becomes a stringent test of the convergence of the calcu- 
lation and thus on the theoretical uncertainties. Further, new potential model calculations 
1T^] and a further comparison with those is worthwhile, given that these new calcu- 



We continue to emphasize 



exist 

lations yield different results at threshold to those in ref. 
the importance of fixing the axial counterterm through a direct experimental measurement. 
This has implications not only for neutrino-deuteron scattering, but also for pp —>■ de^v^,, 
NN — i> NNVu, and parity violating e-d scattering. 



II. THE LAGRANGIAN 

We will briefly review nuclear effective field theory without pions [jTy[. The dynamical 
degrees of freedom are nucleons and non-hadronic external currents. Massive hadronic exci- 
tations such as pions and the delta are treated as non-dynamical, point interactions whose 
effects are encoded in the local operators in the lagrangian. The nucleons are non-relativistic 
but with relativistic corrections built in systematically. Nuclear interaction processes are 
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calculated perturbatively with the small expansion parameter 

^l/a,7,P, |q|) 



Q 



A 



(1) 



which is the ratio of the light to heavy scales. The light scales include the inverse S- 
wave nucleon-nucleon scattering length l/a{< 12 MeV) in the ^5*0 channel, the deuteron 
binding momentum 7(= 45.69 MeV) in the ^5*1 channel, the magnitude of the nucleon 
external momentum p in the two nucleon center-of-mass frame, and the momentum transfer 
to the two nucleon system |q|. The heavy scale A is set by the pion mass m,r- How each 
term in the lagrangian scales as powers of Q can be found in ref. [|19|. It is the nontrivial 
renormalization of the strong interaction operators makes the scaling different from a naive 



derivative expansion |T4|. 



A. The Effective Lagrangian 

The lagrangian of an effective field theory for nucleons can be described via 

C = Ci + C2 + --- , 



(2) 



where contains operators involving n nucleons. Neglecting for the moment the weak- 
interaction couplings, we have 



C, = N^(iD, + ^ 



(3) 



where N is the nucleon field, Mjy is the nucleon mass, Dq and D are covariant derivatives 
and the Dq term is the leading relativistic correction. 

The two-nucleon lagrangian needed for a calculation to NNLO is 



Q ^'\n'^ PiN)\N^ PiN) 



{N^PiNY{N'^(B'^Pi - 215 ■ + Pi^'^)N) + h.c. 



C. 



. AT 

16 

^(^•5o.i)i'/\rT' 



T 



Pi + iD^Pi - 2S'PiEt 



Pi^^ + "D^Pi - 2^0 Pi^ 



N 



C'Q'''°''>{N''PiNy{N^P,N) 



a 



{N^PiN)\N^CD^Pi - 21d ■ Pi^ + P^^^)N) + h.c. 



a 



16 



T 



N \ { N 



P,^2 + S^P. - 2^0 Pi^ 



N 



where Pi and Pj are spin-isospin projectors for the ^Si channel and the ^So channel respec- 
tively. 



P = —(T2(JiT2 , TrP^Pj = -6ij , 

1 

2' 



Pi = —CJ2T2Ti , Tr PjPj = -5ij , 



(5) 



(4) 
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with the T matrices act on isospin indices and a matrices on the spin indices. We incorporate 
isospin symmetry breaking in the C2n^°'^\n = 0, 1, ...) operators, so that C2n^°'^^\ C2n^°'"'^^ 
and C2„'^°'""^ are different. In both the ^5*1 and ^5*0 channels, the strong couphng constants 
C2n have renormahzation scale (/i) dependence. These parameters can be fit to the effective 
range expansion, as detailed in ref. and as reviewed in Appendix ^ 

Relativistic corrections start to contribute to physical quantities at NNLO and have 
generic sizes of 0{p'^/m%) of the leading-order (LO) contribution (see for examples). 



They are suppressed by an additional factor of A^/m^ to other NNLO contributions, and 
thus we can neglect them as small (this is verified numerically). 



B. Weak Interactions 

The effective lagrangians for charged (CC) and neutral current (NC) weak interactions 
are given by 

= -^KJ^ + h.c. , (6) 



= -^IP', , (7) 



where the is the leptonic current and is the hadronic current. We have used Gp = 
1.166 X 10^^ GeV~^. For v-d and V-d scattering, 

ll=vY{l-^,)e , l^, = vr{l-l,)y . (8) 
The hadronic currents can be decomposed into vector and axial-vector contributions 

J; = v^-A- = {y;,-Al)~t{y;-Al) , 

J^^ = -2 sin^ ewV^' + (1 - 2 sin^ OwW^^ - A? - 4 , (9) 

where the superscripts represent isovector components (with S representing isoscalar terms) 
and, later, the currents will be labeled by the number of nucleons involved. 

In a NNLO calculation, the electron mass rrie contributions to the matrix elements are 
counted as higher order (suppressed by a factor of mg/7^), such that q^l^^ = up to NNLO. 
Similarly, if the neutrino mass rriu 7^ but -C 1, then the massless neutrino treatment 

we have here is still applicable. The non-relativistic one-body currents are given by 

A^j!^ = —N^As - gAra)cTkN , (10) 
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where gA = 1.26. Here we have neglected the nucleon vector and axial vector charge radius 
contributions. They only contribute at NNLO with about the same size as the relativistic 
contributions due to the small momentum transfers being considered here. 

We use the notation As for the strange quark contribution to the proton spin 

2S^As = {p |s7^75s| p) (11) 



with the value fSQ] 

As = -0.17 ±0.17, (12) 

where is the covariant spin vector, k^^^ = + k„) and k^-^^ = — t^n) are the 
conventional isoscalar and isovector nucleon magnetic moments in nuclear magnetons, with 

Kp = 2.79285 , Kn = -1.91304 . (13) 

fig is the strange magnetic moment of the proton |^ 



2Mn 

2\ _ Tps/2\ I Tps/2\ 



Up 



Fl{0) = , /x. = G^(0) . (14) 

In ref. ||3^, the SAMPLE experiment found 

^^(-O.lGeV^) = 0.23 ± 0.37 ± 0.15 ± 0.19 n.m. , (15) 

extracted from the proton target experiment. However, the newest deuteron target measure- 
ment suggests that the radiative corrections to the axial form factor were underestimated. 
Thus the central value of G^^(— O.lGeV^) could be 40% smaller Theoretical predictions 
for fig range from -0.8 n.m. to 0.8 n.m.. 

Finally, there are two-body currents relevant to this process. As mentioned in ref. [13], 
there are axial contributions 

= {N^PkNy {N^PaN) - 2^e,JkL2AN^P^N)HN^PjN) + h.c. , 

(16) 

and analogous vector contributions to the deuteron magnetic properties, as described in 



ref. mM 



V^'^ = ee^jL, {N'^P.Ny (V, + ^,) (A^^PaA^) 

-2teL2{N'^P^N)\'V, + ^i){N''PkN) + h.c. . (17) 

The parameters Li and L2 can be fit from np d'y and the deuteron magnetic moment, 
respectively, and are given by Li = 7.24 fm'^ and L2 = —0.149 fm^ at /x = p!9[| . 
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III. u{i?)-D NEUTRAL CURRENT INELASTIC SCATTERING 

For the inelastic scattering process 

u + d^u + n + p, 



(18) 



the differential cross section can be written in terms of leptonic and hadronic tensors /^j^ 



and Wfj,iy as 



Gl |k'| 



■^i(|k'|) F^W, 



^lU 1 



duj'dVt 327r2 |k| 
where ljj{uj') = kQlk'o) represent the initial(final) neutrino energies. 

5i(|k'|) = 1 

for NC processes. The leptonic tensor is given by 

= 8{k''k'" + k^k^" - k ■ k'gf"" + ief^'P^kX) 



(19) 



(20) 



(21) 



The hadronic tensor is the imaginary part of forward matrix element of the time order 
product of two weak current operators. It can be parameterized by six different structure 
functions 



W, 



-Im 



TT 



■/Ae-T(.(P) 



J''(.x)j^(0) d(P) 



P P P^'q^ 



d 



Ml 



M' 



+ iWt 



Ml 



(22) 



where the momentum transfer = k^ — k'^. We can easily see that and don't 

contribute to the differential cross section because qj.^" = 0. In the lab frame (deuteron 
rest frame) the differential cross-section simplifies to 



d'a G\J Ik" 



duj'dil 



27r2 



-Sim) 



9 9 
2Wi sin^ - + W2 cos^ - 



1^3 sin^ 



where the 9 is the angle between k and k and we have used the relation 

9 
2 



q 



-Aujuj' sin^ 



For Vd Vnp scattering, the last terms on the right hand sides of eq.(|2T]) and 
sign. 

Phase space for this reaction is defined by 



Max 



-1,1- 



AMn{v -B)-v 
2uouo' 



21 



< cos^ < 1 



(23) 

(24) 
change 

(25) 
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+ 



+ 



FIG. 1. Diagrams that contribute to Wq and Wi at leading order. The crossed circles denote op- 
erators that create or annihilate two nucleons with the quantum numbers of the deuteron. The dark 
gray blobs arise from scattering corrections involving insertions of Cq operators of the appropriate 
channel, and from the Coulomb Green's function (the light gray blob). The wavy line represents 
the Coulomb interaction, present only in the pp final state. Solid lines represent nucleons, and the 
solid circles represent insertions of electroweak current operators. 



and 



{)<uj' <uj- 2{Mn - \ M% - 72) 



(26) 



where v = u — uj' , Mm is the nucleon mass and B{= 2.2245 MeV) is the deuteron binding 
energy. 

In the deuteron rest frame, 

Wo = Woo = W2-W^ , 
Woi = , 



The structure function can be Q expanded as 

rLO I TxrNLO 



(27) 



(28) 



Now we give the expressions for the structure functions order by order in the perturbative 
expansion. 



A. Leading Order (LO) 

The LO diagrams in fig. |T| give 



LO 



2{C\ 



(1)2 



(29) 
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The vector and axial-vector coupling coefficients are given by 



a 



(0) 



Cf = -As 



sin^ -dw 
1 



a 



(1) 



1 



(30) 



The structure and interaction effects are contained in 

2Mn1 P 



F2 = 



r^2(p,|q|) 



TT (p^ + 7' 

— Im 

TT 



s-^sb, |q|) 



So(p,|q|)'>l?^(p) 



The magnitude of the relative momentum between the final-state nucleons is 2p, with 



P 

where e = 0+. For the NC process, 

S2{P,\CI\) =1 + 

^3(p,|q|) =1- 
Bo{p, |q| ) = ■ 



Mni^ - 'y^ - ^ + ie , 



(31) 



(32) 



q2 (p2 _ ^2^ 
2(p2+^2)2 

q^ (p^ + 37^) 
6 (p^ + 7^)^ 




27r 7 — 



q^ 



12 (7 - tpY 



(33) 



Note that we have further expanded the dependence in powers of q^/ (p^ + 7^) in order 
to later obtain analytic results for the Coulomb contribution in the pp channel. The error 
introduced is numerically small (<^ 1% in total cross section) even though the neglected 
terms are formally LO. 

The NN scattering amphtude has an expansion 



A^A_i + Ao + Ai + --- 
where the subscripts denote the powers in the small expansion parameter Q. 



(34) 



Mn j + ip 
-2Tipd (p^ + 7^) 
(7 + ipY 

-T^ pI + 
(7 + ipf' 



(35) 
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FIG. 2. Diagrams at NLO that contribute to Wq and Wi. The black square represents an 
insertion of the C2 operator in the ^Sq channel, and C2-2, Co,o 'in the ^Si channel. The white 
circle represents an insertion of the two-body currents Li^a and L2,a- Other features are as described 
in fig. ||. 






FIG. 3. Diagrams at NLO that contribute to W^. The light grey circle represents an insertion 
of the nucleon magnetic moment operator. Other features are as described in fig. ||. 



where pd = 1.764 fm. And 



A, 



ip) 



-Att 



M. 



N 



-27r 



1 



'0 p 



M, 



N 



+ ip 



-TT Tq P 



M, 



N 



+ ip 



(36) 



where the scattering length a^^^o,np) _ _23.7 fm is known to high accuracy while the effective 
range Tq '^"'"'^^ = 2.73 fm has a 2% uncertainty [Q. This uncertainty is insignificant in this 
calculation. 



B. Next-to-Leading Order (NLO) 

At NLO, Wq and Wi receive contributions from fig. ^ while receives contributions 
from fig. R|. 
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FIG. 4. Diagrams at NNLO that contribute to Wq and Wi. The black diamond represents an 
insertion of the C4 operator in the ^Sq channel, and C^--^^ C2-1, Co,i in the "^Si channel. Other 
features are as described in fig. |J. 



y^NLO 



(0)2 T-,NLO 



1NLO 



Wi 



NLO 



(1)2 T^NLO 



+ G 



NLO 



w. 



NLO 



o/,^(o)^(o) I P ^ ( riv-)) r^\^) n^'-'n^^'M? \ ^ n^'->) n^'^> T? 

^1<-^A f-^M + + o If-^A f-^M - + "tv^ Oj^j P. 



-2 

,4^(1)^(1) rpLO 



.(i)r^(i) 



(O)r^(O) 



.(l)r^{l)^ 



^(O)r^(O) ipLO 



(37) 



where 



j^NLO 



pLo^j^es^) 

pLO^j^i^So) 



A. 
A, 



{'■Sl)^ 
> 

('^^o)^ 



(38) 



and 



NLO 



(0)2 2MNPd 



Im 



TT 




27 



qNLO 



M 



N 



37r2 



Im 




27r 



i?o(p,|q|) aL?^(p) 
i?o(p, |q| ) fc«ZMAL\^"^(p) + ACfL,,A£',^\p) 



(39) 



Note that Gf'^ = = 0. The L^'s are renormalization scale /^-independent quantities 
defined as 



4:71 (/i - 7) 



L2,A = {p- if L 



2,A 



2nC 



(0), 



(1) 



M 



271 



+ 



Pd 



ip-iY 



(40) 
(41) 



The expressions for C^'^"^ and can be found in Appendix |^. Through W3, we become 

sensitive to weak magnetism at NLO, with coupling coefficients given by 



C'S = -2 sin^ ^^K^^) - Ifis , Ci;^ = (1 - 2 sin^ ^w)^^'^ 



(1) 



(42) 
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FIG. 5. Diagrams at NNLO that contribute to Wq and Wi through two insertions (black squares) 
of the C2 operator in the ^Sq channel, and C2-2-, ^^1,0 ^'^ the ^Si channel. Other features are as 
described in fig. 







FIG. 6. Diagrams at NNLO that contribute to Wq and Wi through mixed insertions of the Li^a 
and L2^A operators (white circle) and C2 operators (black square). Other features are as described 

fig. g. 



C. Next-to-Next-to-Leading Order (NNLO) 



Finally, at NNLO, Wq and Wi receive contributions from figs, 
contributions from fig. |^. 



while W3 receives 



^rNNLO 



Wi 



NNLO 



^NNLO 



I ^r'(O)^ T?NNLO 



NNLO 



4 
3 



NNLO 



_|_ ^(ji^f pNNLO _|_ qNNLO 



2 I -C^^^cf} F^^^ 



I ^^(1)^(1) T^NLO 



(43) 



where 



?NNLO 



^NNLO 



pLO^^Si) 



{'So), 



(44) 



11 





FIG. 7. Diagrams at NNLO that contribute to W3. The filled circle represents an insertion of 
the two-body current Li or L2. Other features are as described in figs. [I| and ^. 



NNLO 



+ a 



V 



Im 



qNNLO 



qNLO 



A, 



Mil 



967r4 

Mn 



Im 



{='■51) .(i5o) 



A, 



('■So)^ 



Im 



7 




(1)' 



67r2 I V 27r 

.(0)7 , ni^)7 \ .(35i) 



|0^?o). 



+4 (cz'L,^A + cr^2,Mj A^r^(p) 

Note that G^^ = G^'"^ = 0. The Lm's are independent quantities defined as 
Stt (/i - 7) 



M 



-4 sin^ i^w (/^ - 7)^ Mw^2 - 27rC 



(1) 



27r ^ 



Af / 



(45) 



(46) 
(47) 



At NNLO, we should mention the effects of other partial waves, beyond 5'-wave. The 
P-wave A^A?^ rescattering does not contribute at NNLO. The D-wave would contribute to W4 
at NNLO, but this structure function does not contribute to the cross-section, so we can 
neglect D-wave initial and final states, also. 

A summary of the expressions in both this and the next section, order-by-order, can be 
found in Appendix O. 



IV. iy{i?)-D CHARGED CURRENT INELASTIC SCATTERING 

For CC processes, a few inputs change from their NC values and there are effects of 
electron/positron mass to consider. For Vd e~^nn, the sign of the last term of eq. ( p3D 
is +. Si, S2, S3, and Bq in eqs. (0), (|23|) and (|3lD still take the same functional forms 
as for NC channels. The phase space is modified, due to the positron mass me and the 
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neutron-proton mass splitting 5m = rUn — nip, to be 
Max 



< cos ^ < 1 



^ 4:Mn{u - B - 6m) ~ u"^ - uj'"^ + m 

— m\ 

me<uj' <uj-2 (^Mn - \J Ml - Mm{B + 5m)^ 



(48) 



In principle, there are also electron mass corrections to eq. (|23|) , but these are only important 
close to threshold and that region will not be probed by SNO. 

For the most part, however, the primary difference between the neutral current and 
charged current cases is the fact that the charged current processes are purely isovector. 
As a result, the charged current results can be obtained from the neutral current structure 
factors with the substitutions: 








c 





\Vud 


V — 




<(1) _ 


Kid 


A — 


V2 



c 



(0) 







C 



(1) 

M 



/9a 
, = 



(1) 







V2Li^A\Vud\ 
1 - 2 sin^ ^ V2\Vud\Li L2 

V — 6m , 



(49) 



where we use IV^^I = 0.975 for this CKM matrix element. The nn scattering amplitude still 
has the same form as eq. (|36D, as do eqs. (|40| ) and (|46|), but with different effective range 



parameters. We have used a^^^o,nn) _ _]^g_5 f^^i and Tq = 2.80 fm. Ref. |3i] indicates 



the uncertainty at the few percent level. It is important to note that a 2% 

uncertainty in a*- •^'o."-") ^^[\\ change the u-d breakup cross section at threshold by 3-4%. 

For Ued — i> e~pp, electromagnetic corrections in the final state are important. But instead 
of solving a three-body problem, we can factor out the Coulomb interaction between the 
electron and two protons by the Sommerfeld factor 



Sim 

Ve 



1 - 
2auj' 



(50) 
(51) 



This approximation is valid because the strength of a single-photon exchange between two 
particles with relative velocity v scales as a/v (Coulomb effect). The effect becomes sizable 
only when the photon exchange becomes nonperturbative, i.e. v a. For an electron, this 
velocity corresponds to a wavelength much longer than the size of the two proton system, 
thus eq. ( ^OD is justified. 

For the proton-proton electromagnetic interaction, the Coulomb contribution is enhanced 
by a factor of 1/f and dominates over other short-distance photon exchange processes. We 
explicitly compute the long-distance Coulomb contribution and encode the short- distance 
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photon effects into the local operators which are then fit to data. In this way we naturally 
incorporate all the isospin symmetry breaking effects in the calculation, except for the un- 
known counterterm Li ^ contribution. Since the Li^a operator only encodes short- distance 
physics, the only contribution to isospin symmetry breaking is through hard photons, and 
thus the effect is of ~ 0{a) and thus negligible. The other symmetry breaking effects on 
Li^A will be estimated later. 

For Ued —>■ e~pp the sign of the last term of eq. (P5|) is negative(— ). The effects of the 



Coulomb interaction are encoded in the S"s and Bq as 



aMN 

V = ^ ■ (52) 



Bo{p, |q| ) = M, 



N 



2% tan' 



d?k \/87r7 ^ii-qk e 
{2nf + 7^ e'^^'T^ - 1 - k"^ + ie 



1 + 



((1 - 2r]l) e - 6T]kk-f - 87 
12(g2 + ^2)2 



2M 



N 



2k 



Finally, for the pp channel, eq. (|36| ) must be replaced by (see Appendix ^ 

-47r 1 



where 



A, 



(i5o,pp) 



(p) 



ip) 



M, 



N 



-2tt 



+ aMNH{ri) 

{i5o,pp) 9 







P 



M 



N 



—IT 



a(^So,pp) 



+ aMMH{ri) 



J'So,PPr A 
'0 P 



M 



N 



ai'So,PP) 



+ aMNH{rj) 



(53) 



(54) 



ZtT] 



(55) 



with ip the logarithmic derivative of the F- function. Note that A^^^^'^'''^ in eq. (|5^ ) is not 
the full pp scattering amplitude. It is the pp scattering amplitude with the pure Coulomb 
phase shift removed, as discussed in Appendix 0. We have used a*- ^O'Pp^ = —7.82 fm and 



,('So,pp) 



2.79 fm. These are known to high accuracy. Furthermore, the values for C^ 
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and 6*2 used in eqs. ( ^OD and (|^) should be replaced by Cq and ^^2'^^^ as discussed 
in Appendix 0. 

Finally, modifications to the phase space and parameters in eqs. (48) and (49) correspond 
to simply changing the sign of 6m, 

5m —>■ ~5m . (56) 

We note again that the expressions from this section, and the previous one, are sum- 
marised in Appendix p. 



V. RESULTS 
A. Unknown Parameters 



As discussed extensively in ref. |T3|, contributions from As, fXs and L2,a which are not 
well constrained are, in fact, negligible (-C 1%) due to quasi-orthogonality between initial 
and final states in the ^Si channel. Thus, up to NNLO the axial two-body counter term 
Li A is the only unknown parameter contributing to each breakup channel. To estimate the 
effect of isospin-symmetry breaking on Li^^, we consider how much Li j^ must vary in order 
for A (defined in eg . (|40|) ) to take on a universal value. This assumes that the symmetry 
breaking effects in Cq , ° , and Li a are all comparable. The effect is ~ 10% in the 
value of Li A at /i = for a natural value of Li a given by dimensional analysis 

I^Ml^^^^5fm^ (57) 

This 10% uncertainty in Li^a corresponds to a 1% uncertainty in the total cross sections. 
This means that we can treat the symmetry breaking effect on Li^a as higher order, and 
that we can take Li^a to be the same in all four channels to the precision of this calculation. 



B. Total Cross-Sections 

We are now able to present a systematic and convergent picture of inelastic neutrino- 
deuteron scattering in all four channels: neutral-current i^x{px)d ~* ^x{px)pf^ (NC) with 
X = e, fi, T, and charged-current — > e~pp and V^d e~^nn (CC). 

We parameterize the cross-sections as 

cr{E) = a{E) + h{E)U^A , 

where the coupling constant of the axial two body current Li^a (with /i = m^) is given 
in units of fm^, and present the results in tables |I| and There are also terms quadratic 
in Li A at NNLO, but they are not significant for values of Li a considered here. We will 
neglect these quadratic terms. 

We have performed this calculation to NNLO largely to test the convergence of the 
calculation and, in turn, be able to place constraints on the theoretical uncertainties in the 
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Testing convergence of the EFT calculation in each of the four channels studied. 
Shown are the relative NLO and NNLO contributions to the total cross-section in each channel, 
for values of Li^a between —5 fm^ and +5 fm^. Note that Vx^px) represents e, /x, orr neutri- 
nos ( antineutrinos ) . 



calculation of i'{V)-d breakup. In fig. ^ we compare the size of NLO and NNLO contributions 
to the cross-sections against the LO contribution. Given the uncertainty in Li^a, we consider 
values —5 fm^ < Li^a < 5 fm^. The shaded areas represent the range of NLO and NNLO 
contributions possible for these values of Li^a- We see that the typical NLO contribution is 
of order 5-20%, while the typical NNLO contribution is less than 5% and, better still, less 
than 3% above 5 MeV. 

A clearer picture of convergence emerges if we decompose the cross-sections into a sym- 
metric piece arising from structure factors receiving contributions at all three orders {Wi 
and W2), and an antisymmetric piece receiving contributions only at NLO and NNLO {W3). 
This is done for the neutral current cross-sections in fig. ^ where we see a cleaner separation 
between size of the NLO and NNLO contributions. From this we can expect, with some 
confidence, that the NNNLO contribution will be less than 3%. This in turn represents the 
formal theoretical uncertainty in our calculation. 



We compare our results to those of the potential model calculations NSGK |12| and 



YHH [|rT| in fig. [T^. We perform a global fit of our results to these model calculations 
with Li^A as the only free parameter. We find that the best fit to NSGK is given by 
LiA^^ = 5.6 fm^, and for YHH L\^^ = 0.94 fm^. These values are both consistent with 
the natural value of Li^a estimated in eq. ( p7D and the quality of the fit is impressive. This 
indicates the 5-10% difference in two potential model calculations is largely due to different 
assumptions made about short distance physics. 

We further examine the ratios of our cross-sections to the potential model results of 
NSGK and YHH in figs. |ll] and [l^ respectively. We see that there are deviations of order 
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E fMeV) 


iy H — ¥ ij 71 Ti 
a f 10-42 cm^l 


( T = P It T} 

b fl0~"42 cm2 /fm^l 


a flO-42 cm2 


U TIT) ( T = P 11 T 1 

) b flO-42 cm2 /fm^l 




00S1 B 


000035 


00312 


000036 


4 


0.0287 


0.00034 


0.0282 


0.00034 


5 


0.0885 


0.0011 


0.0865 


0.0011 


6 


0.188 


0.0024 


0.182 


0.0024 


7 


0.329 


0.0044 


0.318 


0.0043 


8 


0.514 


0.0070 


0.49 


0.0069 


9 


0.744 


0.010 


0.710 


0.010 


10 


1.02 


0.015 


0.968 


0.014 


11 


1.34 


0.019 


1.27 


0.019 


12 


1.72 


0.025 


1.61 


0.025 


13 


2.13 


0.032 


1.99 


0.031 


14 


2.60 


0.039 


2.41 


0.038 


15 


3.12 


0.047 


2.87 


0.046 


16 


3.69 


0.056 


3.37 


0.054 


17 


4.31 


0.066 


3.91 


0.064 


18 


4.97 


0.077 


4.49 


0.074 


19 


5.69 


0.089 


5.11 


0.085 


20 


6.47 


0.102 


5.76 


0.097 



TABLE I. Results for the neutral- current cross-sections, calculated to NNLO. The total 



cross-section for each channel is parameterized as cr{E) = a{E) + b{E)Li^A- 
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E (MeV) 


a (10-42 p^2) 


h (10-42 cm2/fm3) 


a (10-42 cm2) 


b (10-42 cm2/fm3) 


2 


0.00344 


0.000031 






3 


0.0439 


0.00044 






4 


0.146 


0.0016 









U.OZU 


U.UUoD 


U.UZD4 


U.UUUoU 
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0.574 


0.0067 


0.110 


0.0014 


7 


0.914 


0.011 


0.261 


0.0034 


8 


1.34 


0.017 


0.482 


0.0065 


9 


1.87 


0.024 


0.776 


0.011 


10 


2.48 


0.032 


1.14 


0.016 


11 


3.20 


0.042 


1.58 


0.023 


12 


4.01 


0.054 


2.09 


0.032 


13 


4.93 


0.067 


2.68 


0.041 


14 


5.95 


0.082 


3.33 


0.052 


15 


7.08 


0.098 


4.06 


0.065 


16 


8.31 


0.12 


4.85 


0.079 


17 


9.66 


0.14 


5.71 


0.094 


18 


11.12 


0.16 


6.63 


0.11 


19 


12.70 


0.18 


7.63 


0.13 


20 


11.39 


0.21 


8.G8 


0.15 



TABLE II. Results for the charged- current cross-sections calculated to NNLO. The total 



cross-section for each channel is parameterized as cr{E) = a{E) + b{E)Li^A- 
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FIG. 9. A further test of the convergence of the EFT calculation, using the NC channel. The 
left-hand graph arises from comparing contributions from the sum of the vd and Vd channels, or 
specifically the contributions from Wi and W2 which received contributions from all three orders 
of perturbation theory. The right-hand graph looks at the difference between the two channels, or 
specifically the contributions from W3 which appear at NLO and NNLO only. As in fig the 
shaded areas represent the effect of varying Li^a between —5 fm'^ and +5 fm'^. 
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FIG. 10. Inelastic u{jy)d cross- sections as a function of incident vij^) energy. The solid curves 



in the upper graph are NSGK [l^j, while the dashed curves are the EFT results at NNLO, fit with 
5.6 /m^. The solid curves in the lower graph are the results of YHH [lA], and the dashed 



curves are the NNLO EFT results fit with Li^a 
lie right on top of the solid curves. 



0.94 fm'^. In both graphs, the dashed curves all 



5% between our result and YHH in all four channels. The fluctuations seen are in the results 
of YHH, and disagreements at threshold are most likely due to differences in the effective 
range parameters that can be associated with each calculation. However, the agreement 
between our result at NSGK is quite impressive - better than 1% over the whole range 
of neutrino energies studied. For comparison, we note the effective range parameters used 
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FIG. 11. Ratios between the EFT calculation at NNLO and the potential model result of 



NSGK [lij. The upper graph compares the two channels of u-d scattering, and the lower graph 
the two channels of u-d scattering. Agreement is better than 1% over the whole range of energies 
shown, for a single value of Li a = 5.6 frn^ . 







ro^^'P'' (fm) 




To (fm) 




r^^"'"'' (fm) 


Pd (fm) 


This Work 


-7.82 


2.79 


-18.5 


2.80 


-23.7 


2.73 


1.764 


NSGK 


-7.815 


2.78 


-18.5 


2.83 


-23.73 


2.69 


1.767 



TABLE III. Effective range parameters as used in our work and NSGK ^J^/. 



here and those calculated from the potential used in NSGK fSSf (using the Argonne vis 
potential , but with only Coulomb electromagnetic interactions) in table fTT[ 

Of import to SNO is the ratio between charged and neutral current cross-sections 



R 



O'CC 
O'NC 



(58) 



As shown in ref. ||T3|, this ratio R at NLO was insensitive to the value of Li^a, and that is 



still true at NNLO as seen in fig. [T^. We consider two ratios, for both the u-d (relevant to 
SNO) and V-d channels. The latter was the only channel discussed in ref. |T^. Variations 
of Li^A over a large range, from —20 to +40 fm^, leads to a 6% variation in R. This likely 
represents an extreme variation, as it leads to as much as a 90% change in the total cross 
sections. The actual uncertainty in R is almost certainly much less. Further, we can see 
that our values of R agree well with the potential model results (fig. |1^), with the worst 
agreement being seen when comparing to YHH. We use a median value of Li ^ = 3.7 fm^ 
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FIG. 12. in fig. 11, hut comparing to the potential model result of YHH [Hj. The agreement 
is not as good here, using a best-fit value of Li^a = 0.94 frr? . 



for this comparison. The poor agreement with YHH is biased towards threshold, and the 
hkely reasons for this have aheady been discussed. 



VI. SUMMARY AND CONCLUSIONS 

We have performed a calculation to NNLO of all four channels of vd and Vd breakup. 
Our work agrees very well with the latest potential model calculations, subject to the fitting 
of a single unknown counterterm Li^a- Working at NNLO has allowed us to determine that 
our calculation indeed converges at the neutrino energies of interest, which in turn allows us 
to determine a formal theoretical uncertainty of 3% in our calculation. The only outstanding 
issue continues to be the determination of Li a- 

It is imperative that an experimental determination of this counterterm be made. The 
theory without pions cannot reach energies that would allow this to be done with SAMPLE, 
but a breakthrough in the treatment of pions at NNLO might make access possible in the full 
theory. For now, ORLAND offers the best hope, with plans for a high-precision measurement 
in the v-d CC channel. 
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FIG. 13. CC to NC cross-section ratios for v — d and u — d scattering, as functions of in- 
cident energy. Shown are ratios of the EFT NNLO results with Li^a = —20 fm^ (dashed) and 
Li^A = 40 fm^, demonstrating the insensitivity of the ratio to the value of Li^a- 
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APPENDIX A: FITTING THE PARAMETERS OF EFFECTIVE FIELD THEORY 

Here, we summarize the fits to parameters in the theory without pions. 
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FIG. 14. The 'ratio of ratios' between the EFT results at NNLO and those of NSGK [l^j 
and YHH [11], with the EFT results using a median value of Li ^ = 3.7 frr? . The upper graph 



compares the ratios for u — d scattering, and the lower graph the ratios for v — d scattering. Again, 
the agreement with NSGK is excellent. 

1. The ^Si channel 



In the ^Si channel, we use the effective range expansion about the deuteron pole. Thus, 



pcot^o = -7 + ij^Pdiv^ + 7^) + 



(Al) 



To keep the deuteron pole position unchanged at each order in perturbation theory, we 
expand 







+ "-"CO 






— ^2-2 


+ '-'2,-1 


+ ... 




— '-^4,-3 


+ ... 





(A2) 

where the first index continues to denote the momentum dependence while the second in- 
dex indicates the explicit power counting in the Q-expansion. Using the definition of the 
scattering amplitude 



A 



An 



1 



Mjv p cot 5o — ip 



(A3) 
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along with the amphtude computed using the power divergent subtraction (PDS) scheme 
proposed by KSW [|ll 

AeFE = — -j TT (A4) 



withC = EC^L^^y^ 

Matching terms order-by-order in a Q-expansion, one obtains 



'-'0,-1 



-47r 1 



Mn (/i - 7) 
{35i) _ 27r pd72 



(35^) _ 27r prf 
^'-^ " (p - 
^(^51) _ 27r p2^2 



M;v(/x-7)' 



„2 



^^•-^ " (p - 7)3 ' ^^^^ 



where we have neglected relativistic corrections. 



2. The ifi'o channel 

We will deal with the pp channel separately in the next section. For the np and nn 
channels, the procedure is somewhat simpler. The effective range expansion is given by 

pcot(5o = + ^rop2 + ... (A6) 



Here, one obtains 
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3. The ^5*0 pp channel 



In this section we show how we fix the pp scattering parameters to phase shift data 
through matching on to effective range expansion. 

The S-wave pp scattering amphtude can be decomposed into 

A = Ac + Asc , (A8) 

where Ac is the pure Coulomb interaction amphtude with strong interaction "turned off" 
and Asc is the remaining part with both strong and Coulomb interactions. Phase shifts 5 
and 6c are defined by 

dn e^'^ - 1 
~ Mat 2ip 

4:71 e^'^c - 1 , , , 

= m;^^^ • (^^' 

Then 

Mn V 2«p / Mat p (cot Osc — v 
where 5sc ^ 5 — 5c- The Effective Range Expansion states that 



p (cot 5sc -i)=[ + Y^'^p" + v^^'h' + ■■■)- aMMH{ri) 



(All) 



where the shape parameter. This expansion is related to Asc/^^^^'~' which can be 

thought of as the pp scattering amplitude with the pure Coulomb phase shift removed. The 
sum of the diagrams give 



^SC_ ^ ^(i5o,pp) ^ ^ (A12) 



e 



where 



C = Y,ct^\p^-aM^^^r , (A13) 



n=0 



J„ = (,„ + 1 _ _ 1^ _ '^H(„) (A14) 



and 



All \ aM^ 2 y An An: 

with the Euler's constant -^e = 0.577. In eq. ( |A14D , a 4-dimensional pole aM^/(47r(4 — (i)), 
along with the 3-dimensional pole, are subtracted from Jq using the PDS prescription. 

It is convenient to insert the expansion parameter e into eqs. (|A10[ - [AT^ ) to keep track of 
the Q expansion. Then the effective field theory parameters Cgjif can be solved by matched 
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on to the effective range expansion order-by-order in e. The insertion of e is done by the 
t r ansf ormation 



'o 
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(A15) 



The assignment of powers in e reflects the powers in Q scahng of those parameters. For 
example, l/a^^^^ and p scale like Q, aM^ scales like while \H{r]) \ < p/i^aM^) scales like 



1/e. The C'^^^ can be represented in a manner analogous that seen in the ^5*1 channel [19 
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Then from ( |A1U -A12) we obtain 
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The solutions of parameters for nn and np^^^°'^ given in the previous section can be obtained 
by taking a = from the above expressions. 



APPENDIX B: SUMMARY OF SCATTERING CROSS SECTION 

In this appendix, we explicitly list the relevant formulae for the four i'(V)d breakup pro- 
cesses. To make the formulae compact, certain higher order terms have been resummed. The 
Q-expansion expressions to NNLO shown in the text can always be recovered by expanding 
in e to 0{e^). Here e is just a device to keep track of the Q-expansion - its value should be 
set to 1 after the expansion is performed. 
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1. u + d u + n + p and u + d u + n + p 



Differential cross section: 

cPa GW Ik'l 
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27r2 



2W^ sm' ^- + W, cos' - T 2^^^^^W^3 sin^ - 
2 2 2 



with — (+) sign for the u{u)d scattering. 
Phase space: 
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Structure functions (as mentioned above, one can obtain the NNLO results by expanding e 
to C»(e2) ): 
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(/^ - 7)' 

-4sin2i?^(/i-7)'M^L2-27rCl';V^ , (BIO) 

evaluated at ;U = m^, with C7q = —3.56 fm^ and C2 ^"'"'^'^ = 6.55 fm^. Further, 

As = —0.17, and fis = -^2,a = are used in our calculation (though the results are not 
sensitive to these choices). 
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2. v + d^e'^ + n + n 



Differential cross section: 
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evaluated at // = m^, with cj'^°'""^ = -3.49 fm^ and c('^°'"") = 6.46 fml 
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2 ^ 2 



|k'| 

Phase space: 
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4MAr(i/ - 5 + rrin - rrip) - a;^ - a;'^ + ml 



2cuv^ 



<cose<l , 



me < a;' < a; - 2(Mjv - \J M% - Mn{B -mn + nip)) 



Structure functions: 
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^(1) _ \Vud\ ^(1) _ \^ud\ ^(1) _ 1^ I (1) 




+e 
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Im 
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27r 



47r (// --) 
Stt {p - 7) 



+ 



Pd 



(/^ - if 

Pd 
{P - if 



evaluated ai p = m^, with Cq 



-3.77 fm2 and G 



{'So,pp) 



7.50 fml 
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